Log concavity of (i + x) m (i + x k ) 

David Handelman 

Let m and k > 2 be positive integers. We determine necessary and sufficient conditions for the 
polynomial P = (1 + x) m (l + x k ) to be strongly unimodal (frequently known as log concave; see 
[S] for an extensive discussion of the combinatorial significance) . Among the surprises is that the 
criterion, m > k 2 — 3, is also the criterion for P to be merely unimodal (that is, for P of this form, 
unimodality implies strong unimodality). 

Let p = X]T=o aiX% ^ e a polynomial with nonnegative coefficients. Then p is unimodal if the 
distribution of coefficients, (do, cti, . . . , a n ) is unimodal, that is, there exists k with < k < n such 
that do < d\ < • • • < Ofc > dfc+i > • ■ ■ > On (in this definition, if the index j is less than zero 
or exceeds n, then a,j = 0). The polynomial p is strongly unimodal if the function i i— )• dj is log 
concave (which amounts to af > d; + iOi_i) and a^di+2 7^ implies Oi+i 7^ 0. 

In the combinatorics literature, strongly unimodal polynomials are referred to as log concave 
(this can cause confusion when the polynomials are treated as functions, in which case log concave 
has a different meaning). The term strongly unimodal dates back to an undeservedly-neglected 
> 1956 paper of Ibragimov [I], wherein he proves (among other things) that a product of a strongly 

unimodal polynomial with a unimodal one is unimodal, and this property characterizes strong 
unimodality. 

A reparameterization [BH] (known as a tilting in [P]) of the polynomial p is the new polyno- 
mial p\(x) := p(Xx) = ^2,ai\ l x % for some positive real A. It is an easy observation from [BH] that 
p is strongly unimodal if and only if for every reparameterization, p\ is unimodal. 

It is known [H] that if p is any real polynomial with no positive real roots, then there exists N 
such that (1 + x) N p is strongly unimodal. So it seemed of interest to determine the smallest choice 
of N when p has particularly spread out roots (if all roots of p lie in the segment | arg z — tt| < vr/3, 
then p is already strongly unimodal; on the other hand, if some roots have small argument, then 
j the N will likely have to be large). 

psj I We use inner product notation to denote coefficients; thus if p = Y1 a i x ' \ then a, = (p, x l ). 
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THEOREM Let m and k > 2 be positive integers, and set P = (1 + x) m (l + x k ). The 
following are equivalent. 

(a) P is strongly unimodal 

(b) P is unimodal 



( I m + k-3 \ ( m + k — 1 \ . c . , , 

P,x 2 ) < [P,x 2 it m + k is odd 



c 



m-\-k — 2 \ / . 

P, x 2 < p x ~ ~ 11 m + k is even 



(d) m > k 2 - 3. 

Obviously (a) implies (b). Since P is symmetric about (m + k)/2, (b) implies (c). The 
implication (c) (d) is not difficult. The bulk of the work involves showing (d) implies (a). 

For the proof of (c) implies (d) and a small portion of the proof of (d) implies (a), we compute 
(P, x J ) for several values of j near the centre of the distribution. Here P = (1 + x) m (l + x h ), 
and j = (m + k — l)/2, {m + k — 3)/2, (m + k — 5)/2, (m + — 7)/2 if m + is odd, and j = 
(m + fc)/2, (m + /c — 2)/2 if m + A; is even. Write j = (m + k — p)/2, where < p < 5 and is an 
odd integer. 

First, suppose that m + k is odd. Set Q(p) = (P,x j )((m + k - 5)/2)!((m - fe - 5)/2)!/m!. 
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Then 



32 

^ ^ (m + k — 3)(m + k — l)(m — k — 3)(m — k — l)(m — k + 1) 

32 

(m + k- 3)(m + k - l)(m - fc - 3)(m - k - l)(m + k + 1) 

32 / 1 1 \ 

(m + k — 3)(m + A; — l)(m — A: — 3)(m — k — 1) \m — k + 1 m + k + 1 J 
64(m + l) 

~ ((m - 3) 2 - fe 2 )((m - l) 2 - A 2 )((m + l) 2 - k 2 ) ' 

32 

~ (m + fc - 3)(m - fc - 3)(m - fc - l)(m - k + l)(m - fc + 3) + 

32 

(m — A; — 3)(m + A; — 3)(m + A; — l)(m + A; + l)(m + A; + 3) 

32 / 1 1 \ 

~ (m - 3) 2 - k 2 \(m - k - l)(m - A; + l)(m - k + 3) + (m + A; - l)(m + k + l)(m + A; + 3) ) 

64(m + l)(m 2 + 2m + 3(P - 1)) 
~ ((m - 3) 2 - k 2 )((m - l) 2 - A 2 )((m + l) 2 - A; 2 )((m + 3) 2 - A; 2 ) 
32 32 

~ Ilto( m - At + 2/ - 3) + nf=o( m + k + 2l-3) 

32 (nto(™ — A; + 2Z - 3) + Ilto( m + k + 21 ~ 3 )) 
= nf= (("» + 2i- 3) 2 -A; 2 ) 

_ 64((m + l)(m 4 + 4m 3 + (10A: 2 - 14)m 2 + (20A; 2 - 36)m + 5A; 4 - 50A; 2 + 45) 
" nLo((™ + 2Z-3) 2 -A 2 ) 

32(m + A;-5) 32(m — fe — 5) 



0(7) 



nLo( m - fc + 2 ^-3) nLo( m + /;: + 2^-3) 

64(m + l)(m 6 + 6m 5 + (21A; 2 - 41)m 4 + (84A; 2 - 204)m 3 + (35A; 4 - 434A; 2 + 463)m 2 + 

II" „('" • ^ • 2/ 3) 
(70A: 4 - 1036A; 2 + 1350)m - 245/c 4 + 7k 6 + 1813/c 2 - 1575) 



nf =0 (m + A; + 2Z-3) 



The bottom lines for Q(5) and Q(7) were found with Maple. In the latter case, the second factor 
in the numerator had to be broken over two lines (fortunately, it simplifies drastically when we 
substitute m = k 2 — 3). We don't require those two until much later in the paper, but this is a 
convenient time to deal with them. 

We first examine Q(3)/Q(l); by symmetry of the distribution of coefficients of P, if the latter 
is unimodal, then Q{3)/Q{\) < 1. We show that this is the case if and only if m > k 2 — 3 (still 
under the assumption that m + k is odd), and moreover, if m = k 2 — 3 (so that m + k is odd), 
then Q(3) = Q(l). 
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(P^X^ 3 ) ^ Q(3) ^ ( m+ l)( m 2 + 2m + 3(fc 2_ 1)) 

|p )X ^j ~ Q(l) ~ (m+l)((m + 3) 2 -P) 

_ m 2 +2m + 3(fc 2 - 1) 
(m + 3) 2 - A; 2 

/ yyi 3 \ / m -\-k 1 \ 

Expanding this (and assuming, as usual, that m, k > 2), we deduce that I P, x 2 J < I P, x 3 J 
if and only if m > k 2 — 3, with equality if and only if m = k 2 — 3. 

For m + k even, computing ^P, x m+ ? 2 ^ ^ ^P, x 12 ^^ (the ratio of one term away from the 
centre to the central term) is very easy, and we obtain 



P,x 



m + k m — k 



^p^m+h^ 2 V m - k -f 2 in - k - 2 

m 2 + k 2 + 2m 
m 2 + 4m + 4 — A; 2 ' 

It follows immediately that when m + k is even, if and only if m > A; — 2 

with equality only when m = k 2 — 2. 

The two cases (depending on the parity of m+k) yield (c) implies (d). The first case also shows 
that when m = k 2 — 3 (which entails that m + k is odd), the four consecutive middle coefficients 
(that is, positions (m + k ± 3)/2, (m + A; ± l)/2) are equal. 

Next we show that when m = k — 3, IP, x 2 < P, x 2 J , equivalently, Q(5) < Q(3); 

more importantly, the ratio is determined exactly. Note that m + 1 divides the numerator of Q(5); 
this allows a simplification. We note that 



m+fc — 3 



Q(5) _ m 4 + 4m 3 + (10A; 2 - 14)m 2 + (20A; 2 - 36)m + 5k 4 - 50k 2 + 45 
Q(3) ~ (m 2 + 2m + 3(P - l))((m + 5) 2 - A: 2 ) 



m 



4 + 4m 3 + (10A; 2 - 14)m 2 + (20A; 2 - 36)m + 5A; 4 - 50k 2 + 45 



m 4 + 12m 3 + (42 + 2A; 2 )m 2 + (20 + 28A; 2 )m - 3A; 4 + 78A; 2 - 75 ' 



Now assume m = k 2 — 3; the numerator factors (courtesy of Maple) as A: 2 (k — 1) (k — 2) (k + 2) (k + 
l)(k 2 + 7) and the denominator as A; 2 (A; - 1)(A; + 1)(A; 2 + k + 2) (A; 2 - k + 2), so 



Q(5) 



Q(3) 



(fc-2)(fc + 2)(A: 2 + 7) 
m=k2 _ 3 ~ (k 2 + k + 2)(k 2 -k + 2) 

k 4 + 3A; 2 - 28 32 



A; 4 + 3A; 2 + 4 A; 4 + 3P + 4' 



• 1 1 79 f t~, m + k-5 \ / m + k-3 \ 

In particular, when m = A: — 3, IP, x 2 J < IP, x 2 I. 

Both Q(5) and Q(7) have common factors of 64(m + 1); when we factor these out and set 
m = k 2 - 3, the result for Q(7) factors as simply A; 2 (A; 2 - 1)(A; 2 - 4) (A; 6 + 14A; 4 - 63A; 2 - 272). Since 



(m + 7) 2 — k 2 evaluates to (k 2 + 4) 2 — k 2 , we deduce 



Q(7) 



Q(5) 



/c 6 + 14A; 4 - 63k 2 - 272 
m=fc2 _ 3 " ((fe 2 + 4) 2 -P)(P + 7) 

_ fc 6 + 14/c 4 - 63/c 2 - 272 _ 128(/c 2 - 3) 



k 6 + 14/c 4 + 65k 2 + 112 /c 6 + 14fc 4 + 65k 2 + 112 ' 

Now we prove (d) implies (a). It suffices to do this only in the case that m = k 2 — 3 because 
multiplication by 1 + x (or any other strongly unimodal polynomial) preserves strong unimodality. 
Maple has verified strong unimodality of P for 2 < A; < 41. For each positive A, define P = P\ via 

P(x) = (l + Xx) m {l + (\x) k ). 

Assume m = k 2 — 3, where k > 42 is an integer. We show that for all A > 0, Pa is unimodal, 
hence P\ (our original P) is strongly unimodal. Symmetry of Pi implies that it is enough to 
show unimodality of P\ for all < A < 1. The first observation is that the sequence ((P, x J )) is 
increasing for < j < Am/(1 + A) and decreasing for j > Xm/(1 + A) + k. Hence we need only 
consider the integers j in the interval (Am/ (1 + A), Am/ (1 + A) + k). It is trivial to verify that P\ 
is unimodal for < A < k{\ + X)/m, so we may assume A > k(l + A)/m. 
Define (for A fixed), 

p(j) := (P, x J+1 — x J ); then, 

A ^ . = Pti)(m-j) = ^ _ ^ _ ^ + ^ + G+T-J ( A(m - j)(m - j + k) - (m - — k + 1) 



Am-(l + A)j-l+ ^- fc; m J (Xm-(1 + X)j + (1 + X)k-1) 



Next we observe that 



Q+T-fc) _ A j - k + s + 1 



s 

k 



I , and moreover, 

^ m — j — 1 + s 



m—j 
m—j + k 



1 



n!=i(^ -3 + 8-1) \ttl{m -j-l + s) 



U k s =A m - j + s ) 

Thus 

k k k 

Mj) U(m -j + s) = (Am - (1 + X)j - 1) f[(m " J + «) + (Am - (1 + A)j + (1 + A)* - 1) - * + 2); 
11 i 

now setting r = j — Am/(1 + A) with < r < fc, 

= (-l-(l + A)r) 11(1^ ->• + «) +((l + A)(*-r)-l)n( I ^L + r-« + 2) 

-*^(i+*)(*-')-i)n(r^ + !L 4 ±£ )-« 1+ * )r+1) n(rTA- r+ ') 



In particular, p(j) > (that is, the coefficients of P are increasing at j) if and only if 

(1 + AXfc-rQ-l > A Th-r + s 
(1 + A)r + 1 ~ \1 >™ +r -s + 2 



1 l+A 

k 



(1) =A- fe fj 



T+A -r + s 

m s—r—2 

1 l+A A 



n(i + ( *- r l 1+ !^- 2/A ) 

1 V l+A A / 



(We use the forms of both lines 1 and 3 of (1) in what follows.) 

Set M = m/(l + A). The left side of (1) is less than k. The right side of the expression (first 
line) is larger than ((M — k)/(MX + k)) k ; since this is monotone decreasing in A, the minimum 
the expression achieves on the interval (k/M, a] (where a is some number less than 1) occurs when 
A = a. 

We evaluate it, assuming k > 42 and 1 > a > 1/2. 
( M — k \ k _ A *(l +i 



Ma + k 



a k exp—k 




> oT k , , 

2M 2 /k 

> -le" 9 / 2 (1 - 10/fe) > -\.75e- 9 / 2 . 

If we choose any a > 1/2 with a < 1 — (5 + ln k)/k, then it is easy to see that .75e -9 / 2 /a k > k. 
Hence p(j) is negative for < r < k, and thus for any A < 1 — (5 + In k)/k, P\ is decreasing on 
the interval Z n (AM, AM + k); in particular, this implies that P\ is unimodal. So we are reduce 
to considering the case that 1 > A > 1 — (5 + In k)/k. The next, highly intricate, stage is to show 
unimodality for A < 1 — 7/k 2 , and finally there will remain only a few special cases to consider. 

Let G(r) be the logarithm of the left side of (1) and H(r) the logarithm of the right side [first 
line]. From lng = ln((l + \)(k — r) — 1) — ln((l + X)(k — r) + 1), we have 

g(r) = - „, ' + \ , + 1+A 



(l + A)(fc-r)-l (1 + A)r + 1 
k 

Set M = m/(l + A). From H = In h = -k In A + YX (H m + s - r) - ln(M - (s - r - 2)/A)), we 
have 

k 

H'(r) = -f2{—^ + 



s=l 



M + s-r M-(s-r-2)/X 



Now suppose that g(r ) < h(r ), so that G(r ) < H(r ). If G — H is decreasing on [r ,r], 
then G — H is negative on this interval, so that g < h on the interval, hence p(j) < 0, and if we 
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can arrange that this hold for r sufficiently close to fc/2, then (P, x- 7 ) is unimodal. We consider 
when the following inequality holds, 

fc / 1 1 

( ' (k-r-^ir+j^) > Z^ i \M + s-r + M-{s-r-2)/\ 

When this holds, {G — H)' < 0, so that G — H is decreasing. We estimate the right side using the 
first few terms of Euler-Maclaurin summation, 

n pn -i -I 

i=o ■ / ° 2 12 

where |i?o| £ (2^) 3 Jo* l/^l> an( ^ ^ nas cons tant sign on [0, n], then signP = sign/^ 3 - 1 . 

Set M = M - r + 1, so that 1/(M + s - r) = V( M o + «), and set n = fc - 1 and 

/(x) = l/(M + x). Then /' = -l/(M + x) 2 , /" = 2/(M + x) 3 , /'" = -6/(M + x) 4 so that the 
remainder term, R , is negative. Thus 

A 1 y*" 1 dx 1 / 1 J_\ 1 / 1 1 \ 

f^M + s-r-Jv M + x + 2 Vm + A;-1 + M y + 12 V (M + fc - l) 2 + M 2 J 
M + fc~l 1 / 1 1 \ 1 / 1 1 

11 M 2 v M + fc-l M y 12 V (M + fc-l) 2 M 2 

and |i? | < 4/(2vr) 3 (l/M 3 - 1/(M + fc - l) 3 ) < fc/20M 4 . 

The second sum is slightly different: Yli VC^ — (s — r — 2)/A) = X^o _1 M — (s — r — 1)/A; 
set Mi = M + (r + 1)/A and / = l/(Mi — x/A). Then all derivatives of / are nonnegative, so the 
remainder term, Pi, is positive. Thus 

A 1 Z 1 ^ 1 dx 1 / 1 1 / 1 1_\ 

M-^l ~ J M 1 - f + 2 l^Mi - O + Mi J + 12A IjMi - ^±) 2 ~ Mf J + 1 

Mi 1 / 1 1 \ 1 / 1 1 \ 

" Aln Mi-^i + 2 ^Mi-^i + MTj + 12A ^(Mi-^i) 2 ~M!J + U 

and < i?i < /c/20Mi 4 . 

Now we estimate the various terms and add them up. 

M + fc - 1 / fc- l \ _ fc- 1 1 / fc- 1 A 2 l / fc- l \ 3 l / fc- l \ 4 
11 Mo " n V + ^o~y ~^<r~ 2 \Mo~J + 3 V^WcTJ ~4 VlWo" 
k-l 1 /fc - 1 \ 2 1 /fc - 1 N 3 



M 2 V M / 3 V M 



Mi - \ Mi 



AMi 2 V. AMi y 3 V AMi 
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A simple computation reveals that 1/Mq — 1/(M + k — l) 2 < 3(k — 1)/Mq when k > 11 and 
similarly l/(M 1 - (k - 1)/A) 2 - 1/Mf < 3(/c-l)/AM 3 for fc at least 11. We expand l/(M + fc-l) 
and l/(Mi — (fc — 1)/A) by their obvious geometric series expansions, e.g., 

1 1 1 / i-^ + M 2 -r^v+...l 



M + k - 1 M„ l + H M ^ M \ M o / V M o 

1 11 1 / fc-l //c-l\ 2 //c-l x3 



Ml - Mi Mi y AMi V AM i / V AM i 

A simple computation also yields that 1/AMf - 1/Mq < 2(1 - A)/M 2 for k > 11 (the left side may 
be negative; this does no harm). 

Next, we just add everything to obtain an upper bound for —H'(r). We also expand 1/Mq = 
(1/M)(1 - (r - 1)/M) = M" 1 £((r - 1)/M>? and similarly for 1/M 1 . The sum of the cubic and 
higher terms in the expansions of the logarithm is bounded above by k /2M , and that of the 
fourth order and higher terms in all the geometric expansions is bounded by (k — 1)/M 4 . 

rw/ x . ( 1 1 \ k 2 -k / 1 1 
- F(r ^H^ + MTj + ^(AMf"W 

(A;-l) 2 (fc-l) 3 \/l 1 A 1 ( k A n \ n 
+ (^ + L -2-^ + L -^J (M^ + A 2 M?J + M^(y + (fc - 1) J +iil 

2 r(l-l/A)-(l + l/A) r 2 1 - A C" 

" M M 2 " M 3 2AP" /c 4 

^ ± / 2+ (l-A)M + C'M 



M V 2A/c 3 fc 4 

where < C <?? 

Set i? = r + 1/(1 + A) (so 1/(1 + A) < R < k + 1/(1 + A)). The inequality (2) thus holds if 

n „ (l-A)M C'M\ w /c 2 - 3 

it(fc-iZ) (2+ — j +— <M- 



2AA; 3 fc 4 J 1 + A 

We view this as a quadratic in i?. Expanding 

M fc 2 -3 / (l-A)M CM 

(2+ (i_^l + ^M " ^TTa) V 1 " 4Afc 3 " W + 

fc 2 -3 / _ (1-A)M \ 
-2(1 + A)V 4A/c 3 y 

Hence we need only examine the (now normalized) quadratic inequality 
,2 tD , ^ 2 -3 A (l-A)M 



^ -^ + 27l-rA)l 1 -^AA^) > ° ; thati8 ' 



1 1 , (1 - A)(A; 2 -3) A, , 3 3(l-A)(fc 2 -3) 

' 2(1 + A) 4A(1 + A) 2 fc 3 



4(1 + a) v & v fc 2 y A(i + A)yy 2(i + a) 
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We first observe that if A > 1 — (In k + 5)/k (as we may assume by earlier argument), then 
sufficient for the right side (the discriminant) to be negative is that 1 — A > 6/(k 2 (l — a/k) for 
any a < 1. Hence if A < 1 — (6 + l/k)/k 2 , the inequality holds for all values of R in the relevant 
interval, and thus Pa is unimodal. If however, A > 1 — (6 + l/k)/k 2 , we reduce the problem to 
calculating a small number of coefficients. 

Now assume that A = 1 — 7/fc 2 where < 7 < 7 and the inequality in (3) does not hold. The 
right side of (3) becomes 



,2\ 



4-2/ 7 2 \8-4"f/k 2 \ k\ k 2 JJ\(2-"f/k 

3-7/2 1 
< h -• 

4 k 

If 7 > 6 + A/k, this is negative, so the inequality holds for all relevant values of R. For 7 < 6 + 4/fc, 
the inequality is valid providing 



< ^3-7/2 + 4/fc 



So now we investigate what happens when this last inequality fails. We recall that Am/ (1 + A) + r 
must be an integer. Now 

Am (1 -j/k 2 )(k 2 -3) „ 1 
+ r = - U- — - + R- 



1 + A 2-^/k 2 2- j/k 2 

1 00 

= R + (fc 2 _ 4 _ 7 + 3^ £ (7/2A . 2) ; 



fc 2 -4-7/2 + ^ 



i=o 

= " 2 " ' ^ +e, 

where — 37/8A; 4 < e < 7/lO/c . If k is odd, this entails that the fractional part of R be the fractional 
part of 1/2 + 7/4 - / where / ~ (~7 2 + 27)/4fe 2 (to within 7/lOfc 4 ). If k is even, the fractional 
part of R is the fractional part of 7/4 — /. 

In the former case, write R = c + 1/2 + 7/4 — /, and k = 2t + 1 (where c and t are integers). 
Then R — k/2 = c — t + 7/4 — /, and the absolute value of this is supposed to be less than 
y/3 — 7/2 + 4/k/2. If k is instead even, then write R = c + 7/4 — / and k = 2t; then R — k/2 is 
exactly the same as in the odd case. 

Set C = c — t — /, so that C is within a tiny fraction of an integer. The inequality (R — k/2) 2 < 
3/4 - 7/8 + 1/k then translates to (C + 7/4) 2 < 3/4 - 7/8 + 1/k, that is 



7 2 



+ (8C + 2)7 + 16<7 2 - 12 - 16/k < 0. 



If C > l/l, then it is at least as large as a number close to 1, and we see quickly that when k > 10, 
the inequality is impossible. If C is close to zero, then 7 < 2.7 (assuming k > 42). If C is close 
to —1, then 7 < 3.85, and finally if C is close to —2 or any smaller integer, there are no solutions. 
The fact that C must be close to either or —1 means we really only have to calculate a few 
coefficients of the original polynomial. 
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If C is close to 0, then c = t, and we can evaluate Am/(1 + A) + r (an integer): 

(i) Ifk is odd. The nearest integer to i2+(l/2)(fc 2 -4-7/2-2/+2e), where R = k/2+-//A-f+e 
and 7 < 2.7; this is (k 2 + k — 4)/2 (second term to the left of the centre). 

(ii) If k is even. The nearest integer to the same expression, where R = c + — f + e and k = 2t; 
this is the nearest integer to k/2 + (k 2 — 4)/2, exactly as in case (i). 

If C is close to —1, then c = t — 1; when = 2i + 1 is odd, R = t — 1/2 + 7/4 — / + e, so 
Am/(1 + A) + r is 4 - 1/2 + 7/4 - / + (/c 2 - 4 - 7/2 + 2f - 2e)/2, which is just (k 2 + k- 6) /2. The 
result for even A; is the same. 

Thus far, we have shown that if A < 1 — 3.85//c 2 , then Pa is unimodal; and if 1 > A > 
1 — 3.85/A; 2 , then the difference of consecutive coefficients in P\ is decreasing up to j = (k 2 + k — 8)/2 
and for j at least as large as {k 2 + k + 2)/2. 

It remains to show that P\ is unimodal for 1 > A > 1 — 3.85//c 2 (notice that if P\ is unimodal 
for all A < 1, then it is unimodal for A = 1). We have that (P\,x J+ — x 3 ) is decreasing for all j 
in the intervals [0, (k 2 + k - 8)/2] and \(k 2 + k + 2)/2, k 2 + k - 3]. Because of the overlap, it is 
sufficient to show ((P\, x 3 )) is unimodal for Am/ (1 + A) < j < k + Am/ (1 + A). 

Since (Pa, x j+1 — x 3 ) is decreasing on the interval [0, (k 2 + k — 8)/2], it follows that the sequence 

(P\,x 3 )f=+ k ~ 8)/2 is unimodal. Denote (k 2 + k - 8)/2 = (m + k - 5)/2 by s. By the Q(5)/Q(3) 
formula, 

(Pi, X S ) < (Pl,X S+1 ) = (Pi,X S+2 ) = (P!,X S+3 ) = (P l5 X S+4 ) > (P!,X S+5 ). 

Hence the sequence ((Pi, x s+l ))f =0 is strongly unimodal, and therefore the sequence ((Px, x s+l ))f =0 = 
(A s A*(Pi, x s+l )f =0 is unimodal (it is of course strongly unimodal, but this is unnecessary here). 
Moreover, if A < 1, the sequence ((Pa, x s+l ))f =1 (beginning at % = 1 rather than i = 0) is mono- 
tone decreasing. As (P\,x 3+1 — x 3 ) is decreasing for j > s — 5, it follows that the entire tail, 
((P\, x s+l ))i>i is monotone decreasing. 

Let v denote the rightmost mode (integer where the maximum is attained). If v = s, then 
(Pa,x j )* =0 is increasing and thus the entire sequence is unimodal (with mode at s). Suppose 
that v < s. Then unimodality of the entire sequence will result if (Pa,x s ) > (Pa,x s+1 ). Since 
(P A ,x s )/(P A ,x s+1 ) = A- 1 (Pi,x s )/(Pi,x s+1 ), and the second factor is just Q(5)/Q(3) evaluated at 
m = k 2 - 3, we see that if A < 1 - 32/ (fc 4 + 3k 2 + 4), then P A is unimodal. 

This leaves the possibility that A > 1 — 32/ (k 4 + 3 k 2 +4) but v < s. We show this is impossible. 
Since u < s and v is the largest mode, (Pa,x s_1 ) > (P\,x s ). Then the ratio (Pa, x s_1 )/(Pa, x s ) = 
A _1 (Pi,x s_1 )/(Pi,x s ), and again the right factor is Q(7)/Q(5) evaluated at m = k 2 — 3. This 
yields 



Q(7) 



Q(5) 



m=k 2 -3 



> A > 1 - — —5 -, so 

fc 4 + 3/c 2 + 4' 



128(/c 2 -3) , 32 . , , 

1 - fcs + 14^ + 65fc 2 + 112 > 1 " fc 4 + 3fc2 + 4 ; eqmvalently 

1 4(fc 2 -3) 



fc 4 + 3P + 4 fc 6 + 14fe 4 + 65A; 2 + 112' 
A quick check shows that this is false merely if k > 3. This completes the proof of the theorem. • 

Addendum 

Let / be a complex-valued function defined on a neighbourhood of the unit circle, that is at least 
C 2 , and such that f(l) = 1. The variance off is defined as V(f) = f"(l) + f'(l) - f'(l) 2 . Define, 
as in [H2], the class of entire functions, 

S = {entire/ : C -> C,/(l) = 1 if \z\ = 1, then |/(z)| 2 < e -ReV(/)|i-*| 2 J . 
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This class is closed with respect to products and uniform convergence on compact sets. It includes 
all real polynomials all of whose roots lie in the segment | arg z — vr | < ir/3. It is known (and easy to 
prove) [H2] that if p is a real polynomial with p(l) = 1, then there exists N such that ((1 + x)/2) N p 
belongs to £. All known examples of real polynomials in £ with no negative coefficients are strongly 
unimodal, so it will be of interest to check what the optimal value of iV is when p = (1 + x k )/2. 
Another surprise: N must be of order k 4 . 

Set / = {l + x) m {l + x k )/2 m+1 . Then V(/) = mV((l + x)/2) + V((l + x k )/2) = {m + k 2 )/A. 
From |1 — e l9 \ 2 = 4 sin 2 9/2, we have that for z = e %e (it is enough to deal with < 9 < it, since / 
has only real coefficients) to 

\f(z)\ 2 = cos 2m 9/2 cos 2 (k9/2) 

e -V(f)\l-z\ 2 = exp (_( m + fc 2) gin 2 Qjy 

Thus the inequality |/(^)| 2 < e _v ^)l 1_2: ' 2 is equivalent (on multiplying by e m+fe2 and then taking 
rath roots, 

e i+fc 2 /-( cos 2 q/2) . | cos(^/2)| 2 / m < exp((l + k 2 /m) cos 2 9/2). 

Consider the expression obtained by replacing the obnoxious term by 1 and the cos 2 9/2 by a (with 
< a < 1), and for convenience, replace 1 + k 2 /m by s (so 1 < s, and typically, s will be much 
less than 2), 

F(a) :=e sa - ae s . 

On the interval [0,1], F has a unique minimum, when e s = se sa ° , that is, at «o = 1 - lns/s, 
and F(a ) = — (s — 1 — In s)e s /s < 0. In fact, if k 2 /m < 1 (as it will be), then s — 1 — Ins = 
k 2 /m - ln(l + k 2 /m) = (k 2 /m) 2 (\/2 - k 2 /3m + k 4 /Am 2 -...)> 0. If k 2 /m is substantially less 
than 1, F(ao) = —k 4 /2m 2 + e, where e < k 6 /3m 3 is usually good enough in what follows. 

Since F(0) = 1, F has a unique zero on [0, 1), ol\, and F(a) > when a < a±; in particular, 
a necessary condition for / to belong to £ is that cos 2 ir/k < ct\ (since if 9 = 2ir/k, the obnoxious 
term evaluates to 1). We can approximate (for k > 42) ct\ by 1 — 2t + 8t 2 /3+0 (t 3 ) where t = k 2 /m 
will be surprisingly small. 

Now consider p{x) = (1 + x k )/2; then V{p) = k/A, and we consider e -V ( p -" 1-2 ' — |p(2;) 2 |, 
that is, g(9) = e ~ ksin e l 2 - cos 2 \k9 j '2) where z = e %e . It is a triviality (from the additivity of V, 
that is, V(ab) = V(a) + V(6), and that (1 + x)/2 belongs to £ ) to see that if 9 satisfies g{9) > 0, 
then |/(z)| 2 < e~ v (f^ 1 ~ z \ . It is routine to verify that g is positive on an interval of the form 
(0,13 = (3(k), and a quick check shows that if k = 42, then (3 = (3.45101 ± 10" 5 )vr/2A;. Since 
increasing k decreases cos 2 f3(k) (easy), it follows that for all k > 42, g is positive on (0, 3A5n/2k). 
Set a 2 = cos 2 3.457r/4fc. In particular, a sufficient condition for / to belong to £ is that a 2 < oi\. 
It is relatively easy to verify that when we permit larger k, the coefficient 3.45 can be replaced by 
anything less than 4 (for sufficiently large k). 

Doing the numerics, the necessary condition includes m > fe 4 /37r 2 ss .034A; 4 and the sufficient 
condition is m > A^k^/ir 2 « .045/c 4 . (As k increases, the difference between the constants 
becomes arbitrarily small.) In particular, the TV = m in the formulation of this problem must be 
of order k A . 

Why is it so much more difficult for (1 + x) m (l + x k )/2 m+1 to belong to £ than for it to be 
strongly unimodal? I don't know. 
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